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The air dose rate in an environment contaminated with and depends on the amount, 

depth profile and horizontal distribution of these contaminants within the ground. This paper in¬ 
troduces and verifies a tool that models these variables and calculates ambient dose equivalent rates 
at 1 m above the ground. Good correlation is found between predicted dose rates and dose rates 
measured with survey meters in Fukushima Prefecture in areas contaminated with radiocesium from 
the Fukushima Dai-ichi Nuclear Power Plant accident. This finding is insensitive to the choice for 
modelling the activity depth distribution in the ground using activity measurements of collected 
soil layers, or by using exponential and hyperbolic secant fits to the measurement data. Better 
predictions are obtained by modelling the horizontal distribution of radioactive cesium across an 
area if multiple soil samples are available, as opposed to assuming a spatially homogeneous contam¬ 
ination distribution. Reductions seen in air dose rates above flat, undisturbed fields in Fukushima 
Prefecture are consistent with decrement by radioactive decay and downward migration of cesium 
into soil. Analysis of remediation strategies for farmland soils confirmed that topsoil removal and 
interchanging a topsoil layer with a subsoil layer result in similar reductions in the air dose rate. 
These two strategies are more effective than reverse tillage to invert and mix the topsoil. 


I. INTRODUCTION 

Dose reconstruction performed after the accident at 
the Fukushima Dai-ichi Nuclear Power Plant (FDNPP) 
showed that in the main regions affected by the acci¬ 
dent, like the evacuated areas, external exposure to radia¬ 
tion from radionuclides deposited on the ground (ground- 
shine) was the most important pathway contributing to 
effective doses [1, 2]. Since the accident the Japanese 
Government has restricted the sale of contaminated food¬ 
stuffs, and the short-lived tellurium, iodine and xenon 
radioisotopes released (^^^™Te, ^^^Te, 

^^^Xe) have decayed to completion. Therefore, the main 
radiological hazard that persists in the environment is 
exposure to groundshine from radioactive cesium (^^"^Cs 
and ^^”^08). 

Groundshine after a nuclear accident tends to de¬ 
crease due to radioactive decay of short-lived isotopes 
(e.g. ^^^Te, and ^^^Gs) and the penetration 

of fallout radionuclides into soil [3]. The Japan Atomic 
Energy Agency (JAEA) and partner organizations have 
been monitoring the environment in North-East Japan 
since the accident in March 2011 under contract from the 
Japanese Government. In particular, the consortium has 
been measuring radiocesium activity depth distributions 
within soil and monitoring air dose rates at locations of 
flat, undisturbed fields [4]. 
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Understanding the relationship between distributions 
of radioactive cesium within the ground and air dose 
rates is vital for tracking radiocesium migration, predict¬ 
ing future dose rates, and evaluating remediation strate¬ 
gies for reducing dose rates. Previous authors have pub¬ 
lished conversion factors between the concentration of ra¬ 
dionuclides within the ground and various air dose rate 
quantities [5-11]. These conversion factors assume spa¬ 
tially constant radionuclide inventories and depth dis¬ 
tributions. To assist the recovery from the Fukushima 
disaster, Satoh et al. [12] developed a calculation system 
to evaluate air dose rates allowing spatially varying ra¬ 
dionuclide inventories, using a method based on summing 
contributions from radionuclides in different volumes of 
soil. 


In this paper we present and verify a tool to calculate 
ambient dose equivalent rates to high precision for arbi¬ 
trary depth profiles and horizontal distributions of ^^^Cs 
and ^^^Cs fallout within soil. We describe the workings 
of the tool and demonstrate the validity of its predictions 
by comparing against monitoring data of air dose rates 
in Fukushima Prefecture. The tool is applied for un¬ 
derstanding reductions in dose rate seen in North-East 
Japan in terms of migration of radiocesium within soil, 
and for evaluating different soil remediation options for 
contaminated farmlands. 
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II. METHODS 


A. Tool to evaluate air dose rates 

The tool calculates ambient dose equivalent rates at 
1 m above the ground, H*{10) (jiSv/h) [13], hereafter re¬ 
ferred to as air dose rates. The tool consists of conversion 
factors between ^^'^Cs and ^^^Cs activity concentrations 
in different cells and layers of soil, and their contribution 
to the air dose rate. This method allows the dose rate 
to be calculated for any radiocesium depth profile within 
soil and horizontal distribution of the activity, to the limit 
of the precision of the discretization of the ground into 
the different soil volumes. 

The geometry considered is the infinite half-space [3] 
and the land surface is divided into cells by a 12.5 by 
12.5 m mesh (Fig. 1). The tool supports up to 149 by 149 
cells horizontally, which equates to an 1862.5 by 1862.5 m 
area of land. Up to 300 soil layers, each 1 mm thick, 
are modelled below each cell on the mesh. Thus the 
maximum depth of radiocesium contamination is 30 cm. 

The half-space geometry is a model for open and uni¬ 
formly flat land. Therefore any natural or man-made ge¬ 
ographical features that could significantly alter the air 
dose rate, such as buildings, hilly topography or dense 
forests, currently cannot be modelled accurately with the 
tool. The model also does not consider the effects of 
ground roughness on air dose rates. These effects are 
most significant when modelling planes of radionuclides 
on the surface of the ground in the half-space geometry, 
such as in the period immediately after fallout deposi¬ 
tion. However, they are negligible when modelling ra¬ 
dionuclides dispersed within the ground after this initial 
weathering period has completed [14], as is the case in 
this paper. 

The input data for an air dose rate calculation consist 
of and ^^^Cs activity concentrations within each 

discrete volume of soil. The calculation for the air dose 
rate performed by the tool is given by 

H*{10)= ^ (1) 

n,i,j,k 


where Av^n,i,j,k (Bq/m^) is the activity concentration of 
radiocesium in the soil volume, and Cn.ij^k (pSv/h per 
Bq/m^) is the activity to dose conversion factor for that 
soil volume. The index n denotes the ^^“^Cs or ^^^Cs iso¬ 
tope, indices i,j denote the cell position of the soil vol¬ 
ume on the mesh, and k indexes the depth of the volume 
for layers numbered down from the surface. As the cal¬ 
culation in Eq. 1 is a simple sum over all the soil volumes 
in the problem, the run-time of the tool on a standard 
desktop computer is about 10 s. 

The conversion factors for all the soil volumes were cal¬ 
culated using the Particle and Heavy Ion Transport code 
System (PHITS) (version 2.64 - Sato et al. [15]). PHITS 
is a Monte Carlo radiation transport code. The conver¬ 
sion factors represent the dose rate at 1 m above the mid¬ 


dle of the central cell on the mesh per unit activity con¬ 
centration within that volume of soil. The density of soil 
was Ps = 1.6 g/cm^ and air was Pa = 0.0012 g/cm^. The 
soil and air chemical compositions followed Eckerman 
and Ryman [16]. The ^^"‘Cs and ^^"^Cs emission spectra 
were drawn from NuDat2 [17]. Note that the ^^^Cs en¬ 
ergy lines in NuDat2 [17] include the contribution from 
the short-lived daughter product In each case 

the source region was scaled to a vertical line and the 
detectors transformed to planes to maximize the compu¬ 
tational efficiency of the Monte Carlo simulations [18]. 


B. Transforming measured activity depth profiles 
for input into the tool 

The tool cannot accurately simulate scenarios where 
there is significant variation in the soil density horizon¬ 
tally across the simulation region, as a constant soil den¬ 
sity (ps = 1-6 g/cm^) was employed in the PHITS sim¬ 
ulations. However, other constant soil densities (ps 7 ^ 
1.6 g/cm^) or soils with varying density as a function of 
depth (i.e. Ps{z), where z (cm) is the depth in soil from 
the ground surface) can be simulated. The solution is 
to transform the depth coordinate of the source activity 
depth profiles using the mass depth [3] . The mass depth, 
C (g/cm^), is defined as 

az)= [%s{z')dz'. ( 2 ) 

Jo 

A contamination depth profile per unit soil mass mea¬ 
sured in a field survey, A^{z) (Bq/kg), can be recast 
using Eq. 2 into a function of soil mass depth 


Am(C) — 


PsiS')dz' 


(3) 


where Psiz) is the density-depth profile of the soil mea¬ 
sured from field samples. Here tildes are used to distin¬ 
guish the depth coordinate in the field, z (cm), from the 
depth coordinate z applicable for inputting data into the 
tool. 

The activity profile as a function of mass depth, Ani(^), 
can be transformed into a function of z applicable for the 
tool’s constant soil density conversion factors, by reap¬ 
plying Eq. 2: 


A„(z) = A^(C/1.6). (4) 

Finally, an activity concentration depth profile (Av(z)) 
for inputting into the tool can be obtained by multiplying 
Ani(z) by the constant soil density 

Av(z) = 1000 • 1.6 • A^(z). (5) 

The factor of 1000 ensures Av(z) has units of Bq/m^. 
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FIG. 1. The geometry of the simulations, showing the discretization of soil into small volumes with variable activity concen¬ 
trations. 


C. Scenarios considered 

1. Dose rates above flat, undisturbed fields 

Field survey teams have measured depth profiles of 
radioactive cesium in soil at approximately 80 locations 
near to FDNPP since December 2011 [19]. The samples 
were taken from sites at wide, flat areas of land and at 
least 5 m from buildings and trees. Soil samples were col¬ 
lected using a scraper plate. This apparatus was used to 
remove individual soil layers with thickness between 0.5- 
3 cm and with increasing depth from the ground surface 
for radiochemical analysis. Properties analyzed included 
the in situ density (ps(z)) and the ^^^Cs and ^^’’Cs activ¬ 
ity per unit mass (^m(^)) of each soil layer. All survey 
data are published online [ 20 ]. 

We used our calculation tool to predict the air dose 
rate at each sampling site based on the soil activity mea¬ 
surements. We then compared the results with 1 m air 
dose rates measured in the field using hand-held survey 
meters. Data over five soil sampling campaigns were con¬ 
sidered in the analysis. The dates of the campaigns are 
listed in Table I. The measured activity depth profiles 
(Am(5)) were scaled into activity concentration profiles 
applicable for the tool {A^{z)) using the procedure de¬ 
scribed in section IIB. Examples of the processed depth 
profiles for two sites are shown in Fig. 2 (black lines). 

As a scraper plate sample was taken at only one point 
on the ground per location visited in each soil sampling 
campaign, it was assumed that the measured soil activ¬ 
ity depth distribution applied homogeneously across the 
whole region simulated by the tool. Explicitly, the activ¬ 
ity concentrations for each soil layer inputted to the tool 
were identical across all cells on the simulation mesh. 

We also considered a second method for evaluating air 
dose rates from the soil activity samples, based on mod¬ 
elling empirical fits to the activity depth profiles. Mat- 
suda et al. [19] characterized the activity depth profiles as 
a function of mass depth (Am(C)) by fitting exponential 


and hyperbolic secant functions. The exponential depth 
distribution is 

Am(C) = Am,oexp(-C//3), (6) 

where Am,o (Bq/kg) is the activity per unit soil mass at 
the ground surface and /3 (g/cm^) is the relaxation mass 
depth that characterizes the degree of fallout penetration 
into the soil. Figure 2(a) shows a fit of the exponential 
function to soil layer activity measurements at one site. 
The total inventory of contamination per unit area of 
land for this distribution is: 

Amv = 10/3Am,0 ■ (7) 

The factor of 10 ensures Ainv has units of Bq/m^. The 
exponential distribution is a satisfactory model for the 
soil activity depth profile for the first few years after fall¬ 
out deposition [3]. 

Matsuda et al. [19] observed that some of the mea¬ 
sured depth profiles display a maximum in the radioce¬ 
sium concentration below the ground surface. They pro¬ 
posed fitting a hyperbolic secant function to these depth 
profiles, as this function can reproduce a peak in activity 
concentration below the surface. The hyperbolic secant 
function is 

-4m(C) = ^m,o cosh (Co//3)sech(-(( - Co)/^) ■ (8) 

Again (Bq/kg) is the activity per unit soil mass at 
the ground surface, and (g/cm^) is a parameter char¬ 
acterizing the length scale of the distribution. The peak 
in activity occurs at the mass depth Co (g/cm^) below 
the surface. The hyperbolic secant function converges to 
an exponential distribution at large mass depths. Fig¬ 
ure 2 (b) shows a fit of the hyperbolic secant function to 
a measured depth distribution. The total radionuclide 
inventory per unit area for the hyperbolic secant distri¬ 
bution is 

Ainv = 20/3Ai„,ocosh(Co//3)[(7r/4) 

-htan"^ (tanh (Co/(2/3)))]. 


(9) 
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TABLE I. Details of the soil sampling campaigns and sites used for air dose rate predictions. 


Soil sampling Dates Number of Number of sites with 

campaign sites used hyperbolic secant fits 


^St 

Dec 

12-22, 2011 

41 

8 


and 

Apr 17-19, 2012 



2nd 

Aug 

21 to Sep 26, 2012 

82 

13 

3rd 

Nov 

26 to Dec 26, 2012 

81 

28 

4th 

Jun 

4-27, 2013 

80 

12 

5th 

Oct 

28 to Nov 29, 2013 

79 

23 


A^(z) [MBq/m^] 




inputted to the tool for dose rate calculations, (a) Expo¬ 
nential depth profile, and (b) hyperbolic secant depth proHle. 
Site identification codes (IDs) follow Matsuda et al. [19]. 


We followed Matsuda et al. [19] and fitted the expo¬ 
nential and hyperbolic secant distributions to measured 
depth profiles. The hyperbolic secant function was used 
for profiles displaying a peak in activity below the sur¬ 
face, and the exponential function otherwise. Table I 


lists the number of sites in this analysis and the number 
of fits with the hyperbolic secant function. Examples of 
the fitted distributions are shown for two sites in Fig. 2 
(red lines) 

We discounted from the analysis any depth profiles 
showing signs of soil mixing or disturbance [19]. Soil 
disturbance included land cultivation and decontamina¬ 
tion work. We also discounted sites where the air dose 
rate was not measured with a survey meter at the time 
of collecting the soil samples. 

Modelling both the measured and the empirical fits for 
the soil depth profiles yielded two predictions for the air 
dose rate at each site. The air dose rate was calculated 
as the sum of contributions from ^^"^Cs and ^^^Cs, and 
an additional 0.05pSv/h contribution representing the 
background dose rate from natural radionuclides [21]. 

2. Evolution of air dose rates 

In addition to the soil sampling campaigns, JAEA 
and partners have been measuring air dose rates with 
hand-held survey meters at thousands of locations across 
Fukushima Prefecture, including locations with flat, 
undisturbed fields [21, 22]. The monitoring results show 
that dose rates at these locations decreased faster than 
expected by just the physical decay of ^^^Cs and ^^"^Cs [4]. 
Mikami et al. [23] demonstrated that, for the period be¬ 
tween March 2012 and December 2012, relatively little 
migration of the ^^^Cs inventory away from these helds 
occurred. Mikami et al. [21] explained the decrease in 
dose rates, beyond what could be expected by radioactive 
decay alone, by the downward migration of radioactive 
cesium into the soil. 

The relaxation mass depth, /?, characterizes the pene¬ 
tration of fallout into soil for the exponential distribution 
(Eq. 6) . Conversion coefficients published for various val¬ 
ues of /3 can be used to evaluate the Im ambient dose 
equivalent rate given the radionuclide inventory per unit 
area of soil [10]. In contrast, two parameters character¬ 
ize the penetration of the radionuclides within soil for the 
hyperbolic secant distribution - a relaxation mass depth 
P and a mass depth Co for the peak in activity concen¬ 
tration below the surface (Eq. 8). 

To allow direct comparison between exponential and 
hyperbolic secant depth profiles, Matsuda et al. [19] 
proposed an effective relaxation mass depth parameter. 
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(g/cm^), for the hyperbolic secant distribution. The 
effective relaxation mass depth is defined as the value /? 
of an exponential depth distribution yielding the same 
air kerma rate at 1 m as the hyperbolic secant distribu¬ 
tion {K - pGy/h), given an identical inventory of fallout 
radionuclides in both distributions (i.e. Aiav is equal for 
both distributions). 

In this study we used our calculation tool to cal¬ 
culate effective relaxation mass depths for the hyper¬ 
bolic secant fits over the five soil sampling campaigns. 
The effective relaxation mass depths were calculated by 
matching H*{10) from each hyperbolic secant distribu¬ 
tion to an exponential distribution of equal inventory, 

i.e. I3^g (g/cm^). Note that the definition of effective 

relaxation mass depth means that /3^ for an expo¬ 
nential distribution is equal to the relaxation mass depth 
(/3) of the distribution. 

We calculated arithmetic mean, median, minimum and 
maximum /3^g values for each of the five soil sampling 
campaigns (Table II). More sites from the first soil sam¬ 
pling campaign could be used in this analysis than for 
the dose rate predictions (c.f. Table II with Table I), as 
it was not necessary to have a field survey measurement 
of the air dose rate in order to calculate pfg 

We considered the decrement of the components of the 
air dose rate attributable to radioactive cesium fallout, 
i.e. i/*(10) measurements minus a 0.05pSv/h contribu¬ 
tion from natural background radiation, over the first 
four air dose rate surveys [20]. The dates of the air dose 
rate surveys and the mean air dose rates at flat, undis¬ 
turbed fields are listed in Table III. 

We modelled the decrement in dose rates due to ra¬ 
dioactive decay and cesium migration deeper within soil. 

First, we matched the mean values from the soil 

sampling campaigns to the periods of the air dose rate 
surveys (Table III). We then modelled exponential dis¬ 
tributions with P parameters equal to the mean P^g 
values with the tool. The inventories supplied were de¬ 
cay corrected to dates at the middle of each air dose rate 
survey period. The decay corrections assumed an ac¬ 
tivity ratio of released ^^"^Cs and ^^^Cs from FDNPP of 
1.00 on March 11, 2011 [2]. The results show little sen¬ 
sitivity to plausible alternatives (in the range 0.90-1.08) 
for this initial activity ratio. The calculated dose rates 
were then normalized to June 21, 2011, the date at the 
middle of the first air dose rate survey (Table III), for 
comparison with the measured dose rates. As no scraper 
plate soil samples were available for the period of the first 

air dose rate survey (June 4 to July 8, 2011), a P^g 
value of 1.00 g/cm^ was assumed as applicable for this 
period [21]. The sensitivity of the results upon this as¬ 
sumption is checked in the results section. 


3. Spatial variability in soil activity levels 

The calculations with the tool in sections IIC 1 and 
IIC 2 assumed spatially uniform radiocesium distribu¬ 
tions, as only one soil sample was available at each lo¬ 
cation. We considered the effect of spatial variability in 
the radiocesium distribution on evaluating dose rates by 
studying the Ottozawa area. The area lies within 2 km 
of FDNPP and soil samples were taken at multiple loca¬ 
tions across the paddy fields and scrubland in the area. 
Figure 3 shows a map of the area and the soil sampling 
locations. 

This area was remediated between November 2011 and 
May 2012 as part of a decontamination pilot project co¬ 
ordinated by JAEA and is now subject to long-term en¬ 
vironmental monitoring [24]. Remediation consisted of 
removing the top 5 cm of topsoil from paddy fields and 
areas around residential buildings, and cleaning road and 
building surfaces. The air dose rates ranged from 22- 
263pSv/h before decontamination, and dropped to be¬ 
tween 4-110 pSv/h afterwards. Decontamination of this 
area was studied numerically by Hashimoto et al. [25] . 

The soil samples and air dose rates were taken on July 
24, 2014 at 18 locations across the area. Soil samples 
were collected by inserting a cylindrical plastic cup (U- 
8 type, 58 mm internal height, 50 mm internal diameter) 
into the topsoil and collecting the soil contents into plas¬ 
tic bags [26]. One sample was taken at locations 1-16, 
while five samples were collected for locations 17 and 18. 
Locations 17 and 18 lie outside the decontaminated area. 

The ^^^Cs and ^^^Cs depth distributions at locations 
5 and 14 were determined by using a scraper plate to re¬ 
move 1 cm thick soil layers down to a maximum depth of 
10 cm. The activity per unit soil mass in each layer was 
measured using a high resolution gamma spectrometer. 
Unfortunately due to oversight we did not measure the 
in situ densities at the time of collecting the soil sam¬ 
ples. Therefore we had to make an assumption for the 
layer densities. We chose densities equal to the mean 
densities of the soil layers collected over the five scraper 
plate sampling campaigns described in section IIC 1. The 
measured activity depth profiles at Ottozawa were expo¬ 
nential to a reasonable approximation (Fig. 4). 

Scraper plate analyses of the activity depth distribu¬ 
tions were not performed at the other locations (loca¬ 
tions 1-18, excluding 5 and 14). An exponential depth 
distribution was assigned to these locations based on the 
P value applicable at the nearest of locations 5 and 14 
to the site. Locations 1-9 were thus assigned an expo¬ 
nential depth distribution with P = 3.60 g/cm^, and 
locations 10-18 an exponential depth distribution with 
P = 1.83 g/cwP. The total inventory per unit area, 
Ajnv, was inferred by correcting the cylindrical cup activ¬ 
ity measurement for radioactivity at depths greater than 
58 mm as given by the exponential distribution. The in¬ 
ventory for locations 5, 14, 17 and 18, where multiple soil 
samples were taken, was taken to be the mean over the 
various samples. 
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TABLE II. Data for /3^ over the five soil sampling campaigns. 

Soil sampling Number of Number of sites with Mean 
campaign sites used hyperbolic secant fits (g/cm^) 

Median 

(g/cm^) 

Min/i"*‘'°> 

(g/cm^) 

(g/cm^) 

^st 

83 

12 

1.13 

0.93 

0.24 

5.95 

2nd 

82 

13 

1.41 

1.00 

0.11 

8.72 

3rd 

81 

28 

1.56 

1.23 

0.43 

10.36 

4th 

80 

12 

1.64 

1.36 

0.29 

7.73 

5th 

79 

23 

2.17 

1.85 

0.38 

6.41 


TABLE III. Details of the air dose rate surveys and results for the models of dose rate reductions, 
for details. 

‘‘Assumed value - 

Air dose rate survey 

Soil sampling results 

Dose rate measurements 

Models 


Campaign 

Dates 

Campaign 

Mean 

(g/cm^) 

J7*(10) - 0.05 
(pSv/h) 

Relative 

change 

Decay 

only 

Decay & 
migration 

I®' campaign 

Jun 4 to 

Jul 8, 2011 

- 

0 

p 

1.25 

1.0 

1.0 

1.0 

2"‘* campaign 

Dec 13, 2011 to 
May 29, 2012 


1.13 

1.01 

0.81 

0.84 

0.82 

D* part of 3’''^ 
campaign 

Aug 14 to 

Sep 7, 2012 

2nd 

1.41 

0.84 

0.67 

0.76 

0.70 

2"=* part of 3'''* 
campaign 

Nov 5 to 

Dec 7, 2012 

3rd 

1.56 

0.78 

0.62 

0.72 

0.65 

D* part of 4“^ 
campaign 

Jun 3 to 

Jul 4, 2013 

4th 

1.64 

0.64 

0.51 

0.64 

0.57 

2"*^ part of 4**' 
campaign 

Oct 28 to 

Dec 2, 2013 

5th 

2.17 

0.55 

0.44 

0.59 

0.48 




Paddy field 



Truck farm 

^ * 


Wasteland 
Coniferous forest 


|j^' Buildings 


Decontaminated area 


Points of measurement 
Grid for simulation 


Modified from 1:25,000 Topographic Maps ©Geospatial Information Authority of Japan 


FIG. 3. The decontamination boundary and soil sampling locations at Ottozawa. 
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Ay(z) [MBq/m^] 



A/z) [MBq/m^] 



FIG. 4. Measured soil activity depth distributions and expo¬ 
nential fits for the Ottozawa area, (a) Location 5, and (b) 
location 14. 


Two strategies were used to predict the air dose rate at 
locations 1-18. The first strategy assumed that the ra¬ 
diocesium distribution was spatially homogeneous. The 
inventory and depth distribution for that location was 
applied uniformly across the simulation region. 

The second strategy was to model the spatial hetero¬ 
geneity in soil activity levels, as revealed by the soil sam¬ 
ples at the other locations. A 12.5 by 12.5 m mesh was 
overlaid onto a map of the area (Fig. 3). Cells contain¬ 
ing a soil sampling location were assigned the inventory 
and relaxation mass depth for that sample. We adopted 
a simple interpolation method to assign inventories and 
relaxation mass depths to the other cells on the mesh. 
The inventories and /3 values were set equal to the val¬ 
ues applicable at the nearest cell hosting a sampling lo¬ 
cation. Cells equidistant from more than one sampling 
location were assigned inventories randomly from one of 


the equidistant locations. Because of a large disparity be¬ 
tween soil activity levels inside and outside the bounds 
of the remediated area, locations outside the remediated 
area were assigned the inventory of the closest of either 
location 17 or 18. It would also be possible to employ 
other interpolation techniques to assign inventories to 
cells without soil samples, for example based on inverse 
distance weighting techniques or Kriging [27]. 

The assigned inventories for all cells across the area 
are depicted in Fig. 5. The mesh size simulated in the 
tool was 149 by 149 cells for both dose rate prediction 
methods. 


4- Evaluation of farmland soil remediation methods 

To evaluate different methods for remediating farm¬ 
land soils, we used the tool to calculate air dose rates 
after remediation by topsoil removal, reverse tillage, or 
topsoil-subsoil layer interchange. Figure 6 shows a typ¬ 
ical exponential depth distribution for ^^^Cs and ^^^Cs 
within undisturbed farmland soil in Fukushima Prefec¬ 
ture (solid black lines). The ^^'^Cs to ^^^Cs activity ratio 
is applicable on December 01, 2011. This date falls within 
a pilot project on decontamination techniques, and allows 
comparison of dose rate predictions from the tool against 
environmental measurements from the decontamination 
project [24]. 

The relaxation mass depth of the exponential profile 
in Fig. 6 is /3 = 1.13 g/cm^. This follows the result from 
the first soil depth distribution sampling campaign (Ta¬ 
ble II). The air dose rate under these ^^^Cs and ^^^Cs 
inventories and depth profiles is 1.25pSv/h before reme¬ 
diation, including a 0.05 pSv/h contribution from natural 
background radiation. 

The different remediation methods alter the activ¬ 
ity depth distributions of the farmland soil. Figure 6 
shows idealized activity depth distributions after top¬ 
soil removal, reverse tillage, or topsoil-subsoil layer in¬ 
terchange. We used the tool to evaluate the air dose rate 
after completion of each of these remediation options. 

Topsoil removal involves mechanically stripping the 
top 5 cm of the soil, and disposing the excavated soil as 
radioactive waste. The activity profile for the remaining 
soil is, to a first approximation, the exponential distribu¬ 
tion for depths greater than 5 cm prior to decontamina¬ 
tion (dotted red lines in Fig. 6). 

Reverse tillage employs a tractor pulled plough to in¬ 
vert the topsoil. The ploughing creates small ridges and 
furrows on the land surface, which flatten off as the soil 
weathers and relaxes. We approximated the soil as be¬ 
ing homogeneously mixed after this process. Ploughing 
down to a depth of 25 cm thus results in a constant ra¬ 
dioactivity profile initially with depth, followed by the 
exponential distribution at depths below 25 cm (dashed 
blue lines in Fig. 6). 

In topsoil-subsoil layer interchange a layer of topsoil 
is switched with a layer of subsoil. Typically a topsoil 
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inventories assigned to 


cells across the Ottozawa area for simulation of air dose rates at locations 1-18. 
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FIG. 6. Activity depth distributions of ^®'*Gs and ^®^Gs for three farmland soil remediation methods, shown for the in situ 
depth coordinate z. Note break in horizontal axes at 1 kBq/kg to show full distribution of activity with depth. 


layer down to 15 cm is excavated with a digger and this 
soil is placed aside on a plastic sheet. The next 15 cm 
of subsoil is then excavated and stored temporarily on 
adjacent ground. The pit that has been created is re¬ 
filled by first adding a 15 cm layer of the original topsoil, 


and then levelling to the ground surface with the exca¬ 
vated subsoil layer. This strategy can be approximated 
as creating two homogenized layers of activity concentra¬ 
tion below the ground surface. The top layer, down to 
15 cm depth, contains the activity originally between the 
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FIG. 7. Setup of farmland soil remediation simulations: light 
blue area within orange dashed line is remediated land. Land 
in the dark blue area outside is not remediated. Dose rates 
before and after remediation were calculated for the locations 
marked by black spots. 

depths of 15 cm and 30 cm. The subsequent 15 cm thick 
layer below contains the activity that was originally in 
the top 15 cm of soil (green dash-dot lines in Fig. 6). 

To model these remediation scenarios, we considered 
remediation of a 37.5 by 37.5m (1.4km^) area of land, 
equivalent to a 3 by 3 square of cells on the simulation 
mesh (Fig. 7). The simulation models consisted of re¬ 
mediated depth distributions within these cells, while 
the depth distributions outside the area remained un¬ 
changed. Reductions in the dose rates were calculated 
for the center and near to the corners of the remedi¬ 
ated square of land. All dose rate evaluations included a 
0.05pSv/h contribution from natural background radia¬ 
tion. 


III. RESULTS AND DISCUSSION 
A. Dose rates above flat, undisturbed fields 

The predictions for air dose rates above flat, undis¬ 
turbed fields made using the measured activity depth 
profiles compare well with the dose rates measured at 
the sampling sites, as shown by Fig. 8(a). The correla¬ 
tion holds over the range of dose rates covered by the 
dataset (0.09-5.3pSv/h). The predicted dose rates are 
always within a factor of three of the true dose rate, with 
one exception. At one site a dose rate of 0.87pSv/h was 
observed, but the tool predicted 0.075 pSv/h. 

The residual differences between the predictions and 
the measured dose rates are shown in Fig. 8(b). A posi¬ 
tive residual indicates an over-estimation by the tool, and 
a negative residual, an under-estimate. There is no ten¬ 


dency for the tool to either over-estimate or under-predict 
dose rates across the range of dose rates measured in the 
surveys. 

Tyler et al. [28] noted previously that individual soil 
samples can be poor representations of the mean soil ac¬ 
tivity across a wide area. The mean free path in air of 
the primary gamma rays emitted by jg 

around 100 m. Satoh et al. [12] showed that radioactiv¬ 
ity within 500 m contributes significantly to an air dose 
rate. Thus, the total volume of soil contributing to the 
dose rate, down to a depth of 8 cm, is 62 800 m^. As the 
volume of soil collected down to the same depth with 
a 15 cm by 30cm scraper plate is 0.0036m^, the sam¬ 
ple represents only 6 • 10“® parts of the total soil volume 
contributing to the air dose rate. 

Highly variable ^^"‘Cs and ^^"^Cs activity concentrations 
are often found between different soil samples taken at 
the same location. Saito et al. [22] confirmed this was 
the case for soil samples taken in Fukushima Prefecture. 
The variations are caused by heterogeneity in the fallout 
deposition, and by scrubbing and concentration of fall¬ 
out nuclides by local earth surface processes. Therefore, 
a large sampling uncertainty for the inventory of the to¬ 
tal soil volume contributing to the air dose rate should 
be expected if only a single soil sample is available. We 
ascribe the sampling uncertainty from the scraper plate 
measurement as the main source of error in the predic¬ 
tions for the air dose rate shown in Fig. 8(a). 

We next considered the quality of the dose rate pre¬ 
dictions obtained by modelling the empirical fits to the 
measured depth profiles (Fig. 9(a)). The coefficient of 
determination obtained in this case is slightly lower than 
the models employing the measured depth profiles di¬ 
rectly (i?^ = 0.785 versus 0.825). The slight difference in 
values is caused by the predictions for the high dose 
rate locations being slightly less accurate from the models 
employing the empirical fitting functions. The residual 
errors for the predictions at these high dose rate locations 
dominate the squared residuals sum in the calculation of 
R^, and hence the resulting R^ value. 

The residuals for the predictions obtained by modelling 
the empirical fits are shown in Fig. 9(b). Excluding the 
high dose rate locations, the amount of scatter in the 
residuals is comparable to Fig. 9(a). Another way to 
quantify the accuracy of the predictions is to consider 
the mean absolute percentage error. This statistic is less 
susceptible to being skewed by the squared residuals for 
the predictions at the high dose rate locations than R^. 
The mean absolute percentage error of the predictions 
made using the exact depth profiles is 29%. This com¬ 
pares with a mean absolute percentage error of 30 % for 
the predictions obtained by modelling the fitted activity 
depth profile functions. 

The results thus indicate that no significant error is 
introduced by modelling the empirical fits to the activity 
depth profiles instead of the measured step-wise profiles. 
This conclusion necessarily depends on the details of the 
soil sampling procedure. Matsuda et al. [19] measured 
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FIG. 8. Correlation between measured air dose rates from the soil sampling campaigns and simulation predictions using 
measured soil depth profiles as inputs, (a) Measurement-prediction correlation. The dotted line indicates y = x. (b) Scatter 
plot of residual errors in the predictions. The dotted line is y = 0. 
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FIG. 9. As per Fig. 8, except showing dose rate predictions made the exponential and hyperbolic secant fits to the measured 
soil activity depth distributions. 


the activity within 0.5 cm layers of topsoil, followed by 
1 and 3 cm thick layers at deeper depths. If coarser soil 
layer thicknesses are employed, modelling the empirical 
fits may yield more accurate predictions than modelling 
the measured depth profiles, as it is plausible for the 
empirical fits to offer a better representation of the true 
activity profile in the soil. 

B. Evolution of air dose rates 

The distributions of values obtained from the 

exponential and hyperbolic secant fits to the depth pro¬ 


files are shown in Fig. 10(a) for the five soil sampling 

campaigns. Both the mean and median values of 
increase over time (Table II), indicating that the radio¬ 
cesium is migrating deeper into the soil. 

The component of the mean air dose rate attributable 
to radiocesium at the flat, undisturbed fields is plot¬ 
ted for the air survey campaigns with solid diamonds in 
Fig. 10(b). The data are plotted relative to June 2011, 
the date of the first air dose rate survey (Table III). 

The solid line in Fig. 10(b) represents the decrement in 
dose rates that would be expected on the basis of radioac¬ 
tive decay of ^^^Cs and ^^^Cs, and without migration of 
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(a) 




(b) 



Month/Year 


Month/Year 


FIG. 10. (a) Box and whisker plot showing distribution of /3^ values over the five soil sampling campaigns. The whiskers 
show the maxima and minima of the distributions. The boxes show the range in between the 25*^ and TS*** percentiles of 
distributions. The mean values are indicated by asterisks. The full distributions are plotted with symbols, offset to the left of 
each box and whiskers, (b) Measurements and modelling results for the reduction in air dose rate component attributable to 
radioactive cesium at locations of flat, undisturbed fields. The measurements (solid diamonds) show the mean air dose rate 
attributable to radiocesium from the air survey campaigns, normalized to the value at the first air dose rate survey (Table III). 
The vertical bars on the data for the decay and migration model (circles) indicate results when varying /3^ between 
0.5-2.0g/cm^ at the time of the first air dose rate survey (June 2011). 


the radiocesium at the sites. The measured dose rates 
decrease faster than expected by just radioactive decay. 

Mikami et al. [21] explained the additional reduction 
in dose rates between June 2011 and December 2012 by 
migration of the radiocesium fallout deeper into the soil. 
This trend continued through 2013, as shown by the re¬ 
sults of our decay and migration calculations (open cir¬ 
cles, Fig. 10(b)). The decay and migration model re¬ 
sults are reasonably consistent with the measurements, 
although they tend to under-estimate the reduction in 
dose rates by up to 10 %. 

A source of uncertainty in the decay and migration 
model is the choice for mean /3^ for the first air dose 
rate survey campaign (June 4 to July 8, 2011 - Table III). 
The soil sampling campaigns by Matsuda et al. [19] com¬ 
menced in December 2011, so cannot provide measure¬ 
ments to derive a mean value applicable to this 

period. The circles in Fig. 10(b) represent the assump¬ 
tion that /3^ = 1.0 g/cm^ in June 2011. ICRU [3] 

cites j3 values for atmospheric radionuclide fallout in the 
range 0.1-4 g/cm^ for up to one year after fallout de¬ 
position. These results are based on measurements for 
cesium radioisotopes from Chernobyl fallout in Europe 
and Western Russia. 

Takahashi et al. [29] measured depth distributions at 
two grassland sites and three abandoned agricultural 
fields in Fukushima Prefecture between June 21-28, 2011. 
They found that the exponential distribution was a good 
fit for the measured depth profiles, with j3 values in the 


range 0.60-3.08 g/cm^. However, they noted that the site 
giving the highest relaxation mass depth (3.08 g/cm^) 
was pasture land where the soil had been disturbed by 
cattle grazing. Excluding this site from their dataset 
yields a mean /3 value from four sites of 1.20 g/cm^. 

To determine the sensitivity of our decay and migra¬ 
tion model on the choice for the mean /3eff value for 
the first air dose rate survey, we considered the effect 
of varying this parameter in the range 0.50-1.20g/cm^. 
This is a range of values that we consider credible for the 
period between June 4 and July 8, 2011, based on the 

previous literature cited and the mean /3^ value of 
1.13 g/cm^ derived from the first soil sampling campaign 
in December 2011 (Table II). The effect of varying the 

initial value of /3^ in this range is shown by vertical 
bars around circle markers in Fig. 10(b). The ranges in¬ 
dicated by these bars include the measurements, but do 
not permit the decay only explanation for the reduction 
in dose rates. The sensitivity analysis is thus consistent 
with the conclusion that migration of cesium deeper into 
the ground was the main cause behind the additional 
decrement in dose rates. 

There are two other factors that could plausibly ex¬ 
plain the underestimation of the true dose rate reduc¬ 
tion by the model for decay and migration deeper into 
soil (Fig. 10(b)). Although Mikami et al. ]23] suggested 
that little migration of the radiocesium inventory in the 
horizontal direction had occurred, within uncertainties 
their data are consistent with a possible small amount 
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of horizontal migration (on the order of 5-10 % of the 
inventory). 

Another factor is as follows. Although the sites featur¬ 
ing in the air dose rate surveys were chosen to be flat, 
open spaces [21], certain sites may include urban areas 
or areas with roads and paved surfaces at the periph¬ 
ery. The wide field of view of environmental radioactivity 
means that radiocesium within these areas contributes to 
the air dose rate. As the radiocesium within these areas 
has a shorter ecological half-life than areas of glassland 
or agricultural areas [30], i.e. the radiocesium is more 
easily washed away, this could contribute to the under¬ 
estimation of the dose rate reduction by the models in 
Fig. 10(b). 


C. Effect of spatial variability in soil activity levels 
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The Ottozawa area was used to study the effect of spa¬ 
tial variations in the radiocesium distribution on air dose 
rates. The range of dose rates measured at Ottozawa 
in July 2014 varied between 3.5-41.4pSv/h (Table IV). 
This is a higher range of values than measured in the five 
soil sampling campaigns (section Ill A), as Ottozawa lies 
closer to FDNPP than the sites in the five soil sampling 
campaigns and is more highly contaminated with fallout 
from the accident. 

Fig. 11 shows two sets of predictions for the air dose 
rates from soil activity measurements, plotted against the 
dose rates measured in the field. The predictions shown 
in Fig. 11(a) did not account for the spatial variations in 
soil activity levels. Figure 11(b) shows predictions from 
the models incorporating the measured spatial variations 
in the radiocesium inventory. 

It is clear that modelling the spatial variations in the 
contamination distribution yields better predictions for 
the air dose rate. Therefore modelling the spatial varia¬ 
tion is the better strategy if multiple soil samples across 
an area are available to include in the dose rate analysis. 

The coefficient of determination is higher for the pre¬ 
dictions modelling the spatial distribution = 0.753) 
than for the predictions assuming a homogeneous radio¬ 
cesium distribution (i?^ = 0.590). The origin of the dif¬ 
ference in the values is traceable to the model account¬ 
ing for the spatial distribution yielding a better predic¬ 
tion for the highest dose rate site, location 17, with a 
measured dose rate of 41.4pSv/h, than the model as¬ 
suming a homogeneous cesium distribution. 

The mean absolute percentage error for the predictions 
taking into account the spatial heterogeneity of the ac¬ 
tivity is 47 %. This result is higher than the « 30 % mean 
absolute percentage error for the predictions for dose 
rates above flat, undisturbed fields (section III A). This 
difference is also observable by comparing the quality of 
the correlation in Figs. 8(a) and 9(a) with Fig. 11(b). 

There are a number of distinctions between the mod¬ 
elling at Ottozawa and the flat, undisturbed fields that 
contribute to a higher uncertainty for the predictions at 


(b) lOOF 


f 

C/!) 


* 

(D 

3 

C/!) 


; Heterogeneous inventory 

/ 

+/ 

/ 

/ 

+ 


10 V 


1 ^ 


T" 


T" 


1 


+ 

/ 

/+ 


= 0.753 


10 


100 


Measured 77*(10) [gSv/h] 


FIG. 11. Correlation between measured air dose rates and 
predictions from soil activity levels at Ottozawa. (a) As¬ 
sumed a spatially homogeneous ^^^Cs and ^®^Cs inventory, 
(b) Spatially varying inventory informed by all the soil sam¬ 
pling locations. Triangles indicate locations 5 and 14, where 
scraper plate samples yielded the depth distribution. 


Ottozawa. There is a higher degree of measurement un¬ 
certainty for many of the soil samples at Ottozawa than 
for the sites visited in the five soil sampling campaigns, 
as the U-8 cup samples collect smaller volumes of soil 
than the scraper plate. This can be shown by examin¬ 
ing the inventories from the multiple soil samples taken 
at locations 5, 14, 17 and 18 (Fig. 12). There is large 
variation between the inventories between the samples 
at each location. The highest variation is seen for loca¬ 
tion 18, where the largest inventory is four times greater 
than the smallest inventory. 

The mean coefficient of variation for the four Ottozawa 
locations is Cv = 0.41. This is larger than the mean 
Cv = 0.36 observed by Saito et al. [22] for locations within 
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TABLE IV. Results of soil sampling and air dose rate predictions for Ottozawa area on July 24, 2014. Bold indicates soil 
samples taken with scraper plate apparatus. Other samples collected with U-8 cup. Italic indicates values of /3 inferred from 
depth distributions at locations 5 or 14. 


Location Inventory /3 Measured 7L*(10) Prediction based on assumption 

(MBq/m^) (g/cm^) (pSv/h) for Cs distribution (pSv/h) 

^^*Cs ^^'’Cs Homogeneous Heterogeneous 


1 

0.643 

1.89 

3.60 

14.6 

5.3 

9.5 

2 

0.154 

0.453 

3.60 

7.4 

1.3 

5.7 

3 

0.472 

1.43 

3.60 

8.2 

3.9 

9.2 

4 

0.294 

0.868 

3.60 

13.0 

2.4 

5.7 

5 

0.567 

1.67 

3.60 

6.9 

4.3 

8.3 


0.483 

1.47 





6 

2.24 

6.61 

3.60 

12.1 

18.2 

19.9 

7 

0.141 

0.419 

3.60 

6.6 

1.2 

5.6 

8 

0.189 

0.566 

3.60 

7.1 

1.6 

7.8 

9 

1.82 

5.52 

3.60 

8.1 

15.0 

15.2 

10 

0.265 

0.811 

1.83 

7.1 

2.7 

7.7 

11 

0.0342 

0.122 

1.83 

4.1 

0.4 

3.0 

12 

0.258 

0.758 

1.83 

3.5 

2.6 

5.6 

13 

0.305 

0.888 

1.83 

5.4 

3.1 

9.3 

14 

0.0714 

0.210 

1.83 

3.6 

0.7 

4.6 


0.0506 

0.145 





15 

0.0391 

0.114 

1.83 

3.5 

0.4 

6.6 

16 

0.665 

1.95 

1.83 

6.9 

6.7 

12.6 

17 

8.78 

26.1 

1.83 

41.4 

98.0 

89.7 


7.88 

23.7 






12.4 

37.4 






9.95 

29.3 






9.68 

29.0 





18 

9.37 

27.8 

1.83 

27.6 

44.9 

38.8 


2.17 

6.34 






3.95 

11.8 






4.42 

13.3 






2.41 

7.22 






a 100 km radius of FDNPP, which are similar to the sites 
visited in the soil sampling campaigns. Mishra et al. 
[31] independently reported a coefficient of variation of 
0.27 between four samples at another site similar to those 
visited in the Matsuda et al. [19] soil sampling campaigns. 

Another factor contributing to the uncertainty for the 
predictions at Ottozawa include the fact that scraper 
plate samples were only taken at locations 5 and 14. The 
depth distribution at other locations had to be inferred 
from these two measurements. It is notable that some 
of the best dose rate predictions obtained for Ottozawa 
were at locations 5 and 14 (triangles - Fig. 11(b)). 


D. Evaluation of farmland soil remediation 
methods 

We used the tool to evaluate the effectiveness of three 
methods for remediating farmland soils for decreasing 
air dose rates (Table V). We calculated the reduction 
in air dose rate at the center and the corner of a 37.5 
by 37.5 m square area of remediated land, and compared 
with field results from a decontamination pilot project in 
Fukushima Prefecture [24]. Also shown in Table V are 


theoretical limits for the reduction in dose rates, calcu¬ 
lated assuming remediation of all the land surface. 

The performance of the topsoil removal and layer inter¬ 
change methods of remediation are similar. Both meth¬ 
ods yield « 65 % reduction in air dose rates for the square 
area of remediated land. These methods are more effec¬ 
tive than reverse tillage, where the calculations indicated 
a « 50 % reduction in the air dose rate. 

Experience from the decontamination pilot project [24] 
suggests a range of dose rate reductions for topsoil re¬ 
moval and reverse tillage. A number of factors affect 
the percentage reduction in dose rates after land reme¬ 
diation, including the size of the area remediated, the 
homogeneity of the remediation actions, and the magni¬ 
tude of the dose rate before remediation relative to the 
natural background dose rate. The remediation parame¬ 
ters, e.g. the thickness of topsoil removed, or the depth 
of ploughing when performing reverse tillage, may also 
have varied slightly. The general correspondence in Ta¬ 
ble V between the predictions from the tool and observed 
results is encouraging. However, the outlined uncertain¬ 
ties and the criticisms raised by Hardie and McKinley [32] 
about the quantification of errors in the decontamination 
pilot project should be borne in mind when interpreting 
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TABLE V. The percentage reduction in the air dose rate after remediation of farmland soils by three different methods. The 
simulation input data (depth profiles, activity levels, etc.) were applicable on December 01, 2011. Full remediation means that 
all 149 by 149 cells on the simulation mesh were modelled as remediated land. 


Remediation 

Observed Simulation results 



method 

results 

[24] 

Center of 37.5 by 37.5 m 
remediated area 

Corner of 37.5 by 37.5 m 
remediated area 

Full 

remediation 

Topsoil removal 

40-70 % 

73% 

65% 

96% 

Reverse tillage 

30-60 % 

54% 

46% 

71% 

Topsoil-subsoil 
layer interchange 

« 65% 

68% 

60% 

90% 



Cl C2 C3 C4 C5 Cl C2 C3 C4 C5 
Soil Sample 

FIG. 12. Variation in ^^^Cs inventory between soil samples 
at Ottozawa locations 5, 14, 17 and 18. SP denotes a scraper 
plate sample, and Cl, C2, etc. denote the U-8 cup samples. 
Asterisks denote soil samples taken on October 10, 2014 and 
inventory decay corrected to July 24, 2014. Cv denotes coeffi¬ 
cient of variation between the soil sample inventories at each 
location. Cv is the ratio of the sample standard deviation to 
the mean. 


this conclusion. 

One advantage of the layer interchange method over 
topsoil removal is that it does not create waste radioac¬ 
tive soil for disposal. However, the fact that the contam¬ 
inated soil remains at the site after remediation, albeit 
below the ground surface, is tempered by the possible 
availability of the radioactive contaminants for uptake by 
crops or vegetation in future. This point may affect the 


viability of farming these lands after remediation if the 
crops or livestock produced approach food safety limits 
for radioactive cesium content. 


IV. CONCLUSIONS 

The simulation predictions for dose rates at flat, undis¬ 
turbed fields from soil activity depth profiles showed good 
correlation with measurements. Little error was intro¬ 
duced by modelling exponential and hyperbolic secant 
fits to measured activity depth profiles. This conclu¬ 
sion necessarily depends on the experimental parameters 
for measuring activity depth distributions. Soil layers at 
least as fine as collected by Matsuda et al. [19] are rec¬ 
ommended if the data are to be used to evaluate air dose 
rates. Simulations of the Ottozawa area demonstrated 
that modelling spatial variations in contamination levels 
improves the quality of dose rate predictions. This ap¬ 
proach is recommended if multiple soil activity samples 
across an area are available. 

The main uncertainty in air dose rate predictions de¬ 
rived from soil samples is due to the sampling uncer¬ 
tainty for the true soil inventory distribution based on 
the limited volume samples. In situ or mobile gamma 
spectroscopy surveys offer a more comprehensive route to 
assess environmental radiocesium distributions, as they 
are subject to much lower sampling uncertainty [3]. The 
results from these surveys could be used to inform inputs 
for dose rate modelling and improve prediction quality. 

Simulations for the decrement in air dose rates seen at 
undisturbed, flat fields in Fukushima Prefecture for the 
first 20 months following the Fukushima Dai-ichi accident 
were consistent with the hypothesis that radiocesium de¬ 
cay and deeper migration in soil are the main responsible 
factors. Simulations of three farmland soil remediation 
methods for reducing air dose rates demonstrated that 
topsoil removal and layer interchange strategies have sim¬ 
ilar levels of effectiveness, and both methods are more 
effective than reverse tillage. 

Techniques for modelling air dose rates from soil activ¬ 
ity concentrations, such as described in this paper, would 
be effective for evaluating air dose rates in future and for 
planning land remediation works. 
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